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Abstract To account for pulsar frequency glitches, it is necessary to use a neutron 
star crust model allowing not only for neutron superfluidity but also for elastic solidity. 
These features have been treated separarately in previous treatments of crust matter, 
but are combined here in a unified treatment that is based on the use of a Lagrangian 
master functon, so that the coherence the system is ensured by the relevant Noether 
identities. As well as the model obtained directly from the variation principle, the same 
master function can provide other conservative alternatives, allowing in particular for 
the effect of perfect vortex pinning. It is also shown how such models can be generalised 
to allow for dissipative effects, including that of imperfect pinning, meaning vortex drag 
or creep. 



1 Introduction 

Almost immediately after the discovery of pulsars, and their identification as rotating 
neutron stars, it was recognised that, to account for their observed frequency glitches, a 
model of the purely fluid type such as is commonly sufficient in stellar structure theory 
would not suffice, and that an adequate "basic picture" pQ would require the use of a 
more elaborate kind of model allowing for the elasticity of the solid crust. In order to 
do this, a suitable category of elastic solid models was developed in which, as well as 
allowing for the effect of General Relativity (which, in the crust layers, is significant 
but not of overwhelming importance) the main innovation [2] was to use a treatment 
based on elastic perturbations not with respect to a fully relaxed local state (such as 
will usually be available in the context of terrestrial engineering) but with respect to a 
local state that is only conditionally relaxed with respect to perturbations that preserve 
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its density (which in a neutron star may be extremely high compared with the density 
at which the solid would be fully relaxed). 

The motivation for the effort needed to construct global neutron star models [3] 
of such a purely elastic type was however diminished to some extent after it became 
clear that solidity alone could not explain the high rate at which quite large glitches 
were actually observed to occur. To account for this, it was generally recognised to be 
necessary to invoke a mechanism involving angular momentum transfer from a rather 
more rapidly rotating superfluid neutron constituent that can flow without resistance 
through the ionic crust material. To describe such a phenomenon at a macroscopic 
(local but not global) level, various multiconstituent fluid models have been developed 
over the years, both in a Newtonian framework [3] and also (not so much for the minor 
improvement in acuracy that it provides in principle, as because it is for many purposes 
actually easier to work in practise!) in a relativistic framework jSJEJEl- 

Work on the application, at the global stellar structure level, of such multicon- 
stituent fluid models has by now been developed, both in a Newtonian framework jH], 
and in a relativistic framework jU], to such an extent that it now seems worthwhile 
to make the further step of incorporating the effects of solidity, without which the oc- 
curence of the actual glitch phenomenon would not be possible. In order to provide 
what is missing ^U] in a purely fluid descripton of such phenomena, the present article 
shows how a suitably chosen master function A can be used to set up a variational or 
more general Newtonian model of a neutron superfluid flowing within an elastic solid 
background, using a formalism that synthesises the separate descriptions of solidity ^JJ 
and multiconstituent fluidity ^21 E3] that are already available in a form adapted for 
this purpose. An analogous synthesis in a relativistic framework is under preparation 
elsewhere [Tl] . 

2 Newtonian action formulation 

2.1 Dynamic variables and background fields 

Our purpose in this Section is to set up a multipurpose Newtonian variational formalism 
that combines the ones that have recently been developed on one hand for the treatment 
of a multiconstituent fluid ^21 E3] and on the other hand for a non conducting elastic 
solid fJJ • The treatment in the present section does not allow for dissipative effects, but 
is nevertheless designed in such a way as to facilitate their subsequent inclusion in the 
Sectional using the approach that has recently been developed for the multiconstituent 
fluid case [To] . 

Let us start by recapitulating some basic considerations for the generic case, in a 
spacetime with arbitrarily chosen coordinates x^ 1 (fi = 0,1,2,3), of an action density 
scalar, A say, whose role is to act as a Lagrangian master function governing the evolu- 
tion of a set of active or "live" dynamical fields, but whose complete specification also 
involves a set of given - passive or "dead" - background fields. 
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The most general infinitesimal Eulerian (i.e. fixed point) variation of such a La- 
grangian will be decomposable as the sum of two contributions in the form 

8A = 6" A + <5*A (1) 

in which <5^A is the realisable part attributable to a physically possible alteration 
of the configuration of the "live" dynamical fields, while <r A is a virtual part arising 
from mathematically conceivable but (in the context under consideration) physically 
forbidden variations of the "dead" background fields that have been fixed in advance. 
In a typical special relativistic application the only relevant "dead" field might just 
be the Minkowski background metric. What we shall be concerned with here however 
is the more complicated Newtonian case in which, as explained in recent preceding 
work [TI], the necessary set of "dead" fields include the rank-3 spacemetric 7^ and 
the associated time covector (which are restricted to be respectively spacelike and 
timelike in the sense of satisfying t^j^ = 0) together with a (Galilean gauge fixing) 
ether vector e M . As well as these indispensible uniform background fields, the set of 
the "dead" background fields that are fixed in advance will be taken here, in the first 
instance, to include the generically non uniform gravitational potential 0, so that the 
complete background variation contribution will be given by an expression of the form 

d a dA c au <9A r OA r „ dA r 1 

For more general purposes, the scalar would need to be promoted from being a given 
background to the status of a "live" dynamical field, which means that the last term 
in (J2J) would be transfered to the "live" variation contribution S^A, which in that case 
[TB*] would also include a term allowing for the dependence of A not just on <p but also 
on its gradient with components <p )fl . 

It is to be remarked that the construction of many of the relevant quantities will 
involve the rank 3 Euclidean space metric that is induced by the choice of the ether 
gauge vector e M acording to the specifications 

r)» v e v = , v„Y P = e% = rf . (3) 



In the more specialised work on the fluid case [13] this uniform space metric was 
denoted by 7^, but the latter will be used here for the variable space metric that is 
to be defined below. As noted in our preceding work [2], the induced variation of the 
uniform space metric @3J) will be expressible (using round brackets to indicate index 
symmetrisation) by the formula 

= -VwVv* h pa - 2tQMv) P , ( 4 ) 

in which the last term shows the effect just of an infinitesimal Galilean transformation, 
under which 8^i pu would simply vanish. 

The particular kind ^2] of background variation in which we shall be interested 
is one that is simply generated by the action of a displacement field say, so that 
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the ensuing field variations will be given just by the negatives of the corresponding Lie 
derivatives. For the background fields that are uniform (in the sense of having vanishing 
covariant derivatives) the resulting variations will be given just by 



fry"" 



5e v = e?V£ v , 



(5) 



while for the gravitational potential we shall simply have 



5(f) 



(6) 



In so far as the "live" contribution is concerned, a further subdivision arises in the 
cases such as those of interest here, which are characterised by a constrained variation 
principle, meaning one whereby the "on shell" evolution condition is that of invariance 
of the relevant action integral with respect to a compactly supported perturbation of 
the dynamical fields that is not entirely arbitrary but constrained by an appropriate 
admissibility condition. In the applications under consideration, the relevant admissi- 
bility condition will be interpretable as the requirement for the change to represent a 
"natural" variation of the same given physical system, whereas more general changes 
would represent a replacement of the system by a slightly different system within the 
same category. In such a case the "live" field contribution will be decomposible (though 
not necessarily in a unique manner) as a sum of the form 



in which #"A denotes a part the would be inadmissible for the purpose of application of 
the variational principle, whereas frA is a "natural" variation that would be allowed for 
this purpose. This means that for any (unperturbed) configuration that is "on shell" - 
in the sense of satisfying the dynamical evolution equations provided by the constrained 
variation principle - a generic admissible variation must satisfy condition 



using the symbol = to indicate equivalence modulo the addition of a divergence (which, 
by Green's theorem, will give no contribution to the action integral from a variation 
that is compactly supported). 

In the kind of medium with which we are concerned, the constrained variables will 
be a set of current 4- vectors n M with "chemical" index label x (which, in the kind of 
astrophysical aplications we have in mind, might include the separate neutronic and 
protonic contributions to the conserved total baryonic current) while the unconstrained 
fields will consist of a triplet of scalar fields q A (a =1,2,3) that are interpretable as 
local coordinates on a 3-dimensional material base manifold, and from which, together 
with their gradient components q A ^, an associated particle number (which might count 
ionic nuclei or crystaline solid lattice points in the relevant microscopic substructure) 
density current, n" say, can be constructed in the manner prescribed below. The 



fr?A = 5^ A + fr*A 



(7) 



5^ = 0, 
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corresponding "live" part of the generic action variation will thus be expressible in the 
form 

6"A = E^Sn* + K S t» + > ( 9 ) 

in which, for each separate current, the partial derivative 7r x = dA/dn£ is interpretable 
as the corresponding 4-momentum per particle. 

It is to be observed that the partial derivative V v = dA/dq A u has the notewor- 
thy property of being tensorial not just with respect to transformations of the space 
time coordinates x u but also with respect to transformations q A —>■ q A of the material 
coordinates, q A whose effect will simply be given by 

This good behaviour contrasts with the comportment of the remaining partial deriva- 
tive, dA/dq A , for which the corresponding transformation has the non-tensorial form 

dA OA dq* mv d 2 q c , . 

h-> — + V q B — . (11) 

QqA QqB Qq\ c QqBQq^ 

For each separate current n M , the admissible variations (as generated [T21 by world- 
line displacements) will be specifiable by a corresponding displacement vector field, £^ 
say, by the formula 

A< = -4£<-<V^ (12) 
in which the Lie derivative is just the commutator 

4£< = ev,<-<v,e. (is) 



The variations of the material base coordinate fields q A will similarly be expressible 
in terms of their own displacement vector field £ M , by an expression of the corresponding 
form 

SV = -q^»i (14) 

but in this case it is evident that no restriction is entailed, so that without loss of 
generality we can take 5^q A = 0, which means that the inadmissible part (if any) of a 
"live" variation will reduce simply to 

5 s A = E x 7r x l/ ^. (15) 

It can be seen (using integration by parts) that, in terms of the displacement fields 
introduced by (J 12)1 and (|14jl . the generic admissible variation for the current carrying 
medium will be expressible -modulo a variationally irrelevant divergence term - in the 
form 

^A--/S,r-£ x /*£x, (16) 
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in which, for each value of x, the covectorial coefficient will be interpretable as the 
non-gravitational force density acting on the corresponding separate constituent, while 
the extra coefficient f s v is the supplementary force density due, as discussed below, to 
stratification or solid elasticity, that acts on the underlying ionic (crystalline or glass 
like) lattice structure. 

If we postulate that the system should obey the constrained variational principle 
that is expressible as the imposition of the "on shell" condition (jHJ), then it is evident 
from (pUjj) that the ensuing dynamical field equations will be expressible simply as the 
separate vanishing of each one of these force densities, i.e. as the requirement that 
we should have f u = and (for each value of x) = 0. We shall however be 
mainly concerned with other possibilities, in particular with cases in which some of the 
constituents, labelled x=c and x= 0, are "confined" in the sense of being subject to a 
convection condition of the form 

<??„ = 0, n fl Y„ = 0, (17) 

which means that they have to move with the underlying ionic flow. When the appli- 
cation of the variation principle is subject to the corresponding convective constraints 

K = tf = e, (is) 

the ensuing system of dynamical field equations will no longer entail the separate 
vanishing of f s u , f\ and f c v , but only of the amalgamated ionic force density f \ that is 
defined as their sum, 

f\ = fl + fl + fl (19) 
which is interpretable as the net force density acting on the integrated system consisting 
of the ions in combination with the convected constituent. 



2.2 Noether identity for Newtonian covariance 

As in the preceding work [12], let us now consider the Noether type identity that is 
obtained (regardless of whether or not the variational field equations are satisfied) by 
taking the variations of all the relevant fields, "live" as well as "dead" to be given by 
the action of the same displacement field whose effect on the background fields was 
given by (jSJ) and ([S|l. so that the resulting effect on the Lagrangian scalar will be given 
just by the corresponding displacement variation 

SA = -f "V„A . (20) 

For the unconstrained fields q A the effect of such a displacement will be given, according 
to (JHJ), simply by 

5q A = -££q A = S^q A . (21) 

On the other hand, for the constrained variables n v ,it can be seen by setting £ u to 
in (|12|) that we shall have 

bnl = -g£n» = + 5«< , (22) 
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in which the extra "unnatural" (variationally inadmissible) contribution will be given 
by 

^K = K^- ( 23 ) 
It can thereby be seen from JIJ) using (jl6j) and ()20)1 that, as an identity which will 
hold, modulo a divergence, for any dynamical field configuration whether or not it is 
"on shell" , we shall have a relation of the form 

C (fl + £ x /* - V„A) = 5 s A + 5* A . (24) 

We now evaluate the terms on the right by substitution of (|2^|) in (|T5|l and by sub- 
stitution of the formulae JHJ) and (JBJ) in (J2J) for an arbitrary displacement field The 
identity (|24j) can thereby be rewritten more explicitly in the form 

e {fl + s x /- - P v,0) = -nv,r - rv^n , (25) 

in which the effective (passive) gravitational mass density can be read out simply as 

<9A , . 

" = -w (26) 

while the corresponding (gravitationally passive) geometric stress momentum energy 
density tensor can be read out as 

r)A BA BA 

n = (a - e/x) K - *r-^ + ~ e "d? ■ (27) 

Since the arbitrary field can be taken to be non zero only in the immediate neigh- 
bourhood of any chosen point, the Noether identity (|23j) implies that, at each point, 
this stress energy tensor will satisfy a divergence identity of the simple form 

V M n = / s , + E x /»-pV^. (28) 

In the "on shell" case for which the variational evolution equations are satisfied, the 
non gravitational force density contributions on the right of (|28|) will drop out, leaving 
just the final term on the right which represents the gravitational force density which 
will always be present since the gravitational potential <fi is treated as a background 
field. 

The tensor T^ v can be employed in the usual way for construction of the energy flux 
vector U u and the energy denstity U according to the prescriptions 

= -T» v e v , U = , (29) 

while the corresponding space momentum flux JJ^ and space momentum density iJ^ 
will be given by the prescriptions 

= T^ au n y = t^n^ , (30) 
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which can then be used to determine a pressure tensor P^ u that is specified by the 
decomposition 

Hi* = P »v + e »n v . (31) 

( It is to be remarked, as a minor caveat, that the systematic notation scheme previously 
developed for cases in which all the constituents were of purely fluid type [T3] deviated 
slightly from the - even more systematic - scheme used here, not only by using 7^ for 
what is denoted here by 77^, but also by using P^ v just as an abbreviation for T£ ■) 
It can be seen from the Noetherian construction (|29|) that the pressure tensor that is 
specified in this way will automatically be symmetric and strictly spacelike, 

pi*> = P ^ ; pp, tv = ? (32) 
since it will be given explicitly by 

8A 

f»» = (a - e x «) y u - z-r-^-r ■ ( 33 ) 

It can also be seen that the 3-momentum density will contain no contribution from any 
part of the Lagrangian that is independent of the Galilean frame, as is the case for the 
purely internal action contribution A int that will be discussed below: we shall simply 
have 

n u = -—Y u => n y =0. (34) 

Q e a I int V ) 

The corresponding formula for the scalar energy density is 

<9A OA 

c , =Exn< _ A _ t _ +e ._. (35) 



2.3 Canonical formulation 

If, instead of working it out in terms of force densities modulo a divergence as in (|2~%Jl 
and (|25|). we insert the displacement contributions (|21j) and (|2*2*|) to the "live" variation 
directly in we get an identity in which the coefficients of the locally adjustable 
fields £ M and must vanish separately. The former of these conditions reduces to 
a triviality, but the latter provides a relation showing that the geometrically defined 
stress energy tensor (|28J) can be rewritten in the equivalent canonical form 

T\ = A5t + S x (ttV^ - 7r-n| n) - V\ , (36) 

in which the stress contribution at the end is given in terms of the coefficient introduced 
in ([TOD by 

V\ = V\q%. (37) 
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The "live" variation formula Q also provides what is needed for the derivation 
of the corresponding expressions for the force densities introduced in (|16|) which for 
the fluid constituents will have the form that is already familiar from the preceding 
work in terms of the corresponding generalised vorticity 2- forms defined, using 

square brackets for index antisymmetrisation, by 

w%, = 2V [M 7r*] , (38) 

these non-gravitational current force densities will be given by 

f^n^ + nV^. ( 39 ) 

while the non-gravitational force density acting on the underlying atomic structure of 
the medium will be given by the prescription 

in which the Eulerian derivative is given by a prescription of the usual form 

— — - VP" V v - — 

5q^~dq^ " a V a-q^T- \ ) 

It is important to notice that the non-tensorial base coordinate transformation property 
(II lj) of the first term in this formula will be cancelled by that of the second, so that 
the effect on the Eulerian derivative of a change of the material base coordinates will 
be given simply by 

5 A 5 A dq B 

q A ^q A -=— i-> -=— -£r- , 42 

Sq^ SqB dq^ 

which shows that the (stratification and solid elasticity) force density (j4*Uj) is invariant 
with respect to such transformations, an important property that is not so obvious from 
its detailed expression 

dA 

fl = Qjzt* + - V„ (7>XJ . ( 43 ) 

in which only the final term is separately invariant. 

It is important to observe that, regardless of what - variational or other - dynamical 
field equations may be imposed, the force density P3*j) will never be able to do any work 
on the medium since, as an obvious consequence of (|4U|) it must automatically satisfy 
the further identity 

fX = , (44) 
where is the unit 4-velocity of the medium as specified by the defining conditions 

q\u v = , t v u v = 1 . (45) 
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If the variation principle JHJ) is imposed, so that - as remarked above - by ()16|) the 
separate force densities f „ and will all have to vanish, then it evidently follows from 
the Noether identity (|28|) that the divergence condition 

V M n = -pV^, (46) 

must hold. This condition is interpretable as an energy-momentum balance equation 
that must be satisfied whenever the system is not subject to any non-gravitational 
external forces. So long as we are concerned only with a system that is isolated in this 
sense, so that the energy-momentum balance condition (|4*U|) is satisfied, it follows from 
that, even for a more general model admitting the presence of dissipative or other 
internal forces, these must be be specified in such a way as to satisfy the total force 
balance condition 

f v + E x fl = 0, (47) 

in order for the model to be self consistent. More particularly, in view of f)44j) it can 
be seen that the fluid force densities by themselves will have to satisfy the condition 

U "E x /*=0, (48) 

expressing energy balance in the local rest frame of the underlying medium. 



3 Application to superfluid neutron conduction in 
convective solid 

3.1 Implementation of the confinement constraint 

In the appendix of the third part ^3] of the preceding series of articles on multicon- 
stituent fluid models in a Newtonian framework, an account was provided of the kind of 
application that provided the main motivation for this work, namely a model in which 
the chemical variable x runs just three values, which are x=f and x= c say, labelling 
a "confined" baryon current and a "free" baryon current n^, together with a third 
value, x=0 say, labelling an entropy current s v = suX, which (because of the relatively 
low temperatures involved) will play a relatively minor role in the pulsar applications 
we have in mind. The new feature in the present treatment is allowance for the solidity 
that (again because of the relatively low temperatures involved) will characterise the 
"confined" constituent, whose flow will therefore be constrained to be aligned with the 
flow vector u u of the solid structure as introduced above, so that it will have the form 

n v c =n c u v . (49) 

As in the previous treatment ^3] we shall restrict our attention here to the case in 
which the (in any case relatively unimportant) entropy current is convected with the 
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solid structure, so that the thermal rest frame will be identifiable with that of the solid, 
Uq = u v , which means that the entropy current will be constrained to have a form 
analogue to namely 

s u = su v . (50) 

The part of the baryon current that is confined according to ()4Dj1 (meaning that has 
to be comoving with the underlying atomic or crystal structure of the medium) will 
add up with the "free" part to give a conserved total baryon current 

< = < + <, V„< = 0- (51) 



It can be seen to follow from this current conservation requirement that, in terms 
of the "free" current vorticity 2-form w v , the energy conservation identity (|48|) will 
reduce to the form 

u»w\ v n v i = (£ c - £ f ) V„< - QW v s v , (52) 

in which the relevant rest-frame particle energies, £ c and £ f , and the temperature, 
say, are given in terms of the corresponding 4-momenta, as introduced in © , 
by the definitions 

S C = , £ f = 7T f „ , = k\ . (53) 



It is to be noticed that these definitions ()53|) are not affected by the freedom to adjust 
the specifications of the confined particle 4-momentum covector tt c and the thermal 
4-momentum covector ift due to the constraints (J4*9|) and (|50jl on the corresponding 
currents n" and s", of which the former can be seen to have the form (|17|) which 
evidently means that the specification of q\ in (jHJ) will also be ambiguous. These 
ambiguities will however need to be resolved in order for us to be able to proceed in 
a well defined manner. To obtain a formalism that matches smoothly with what has 
already been set up for the fluid limit ^3], we need to remove the ambiguity in 7r c „ by 
the imposition on 7r® and q A v of appropriate conditions, which can be taken to be that 
the former should only have no space component, which means that will be given in 
terms of the temperature by 

< = -et #4 , (54) 

while the latter should have only space components, meaning that it will satisfy the 
condition 

V\ t v = . (55) 
The stress energy tensor (J36)) will thereby be obtained in the form 

7* = (A + Bs - vr f p < - ttX) *"„ - + + « - V» A q% . (56) 
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(Instead of the convention (J55j) . an obvious alternative ansatz would be to use the 
analogue of (p>lj) to the effect that n c should have no space component, so that it 
would take the form 7r c „ = —£ c t tl . Such an alternative choice would be adequate for 
the specialised purpose of matching the formalism used for the non-conducting solid 
limit case [TT], but not for broader objective of consistency with the formalism developed 
for the multiconstituent fluid limit case fT2\ FEU ITS].) 

It is computationally convenient and physically natural to take the material base 
space to be endowed with a measure form say that is specifiable in terms of antisym- 
metric components, ^i ABC say, that are fixed in the sense of depending only on the q A , 
and that will determine a corresponding scalar space-time field n\ by the determinant 
formula 

ni = t^Ubc^Idef 7 ad 7 se 7 cf , (57) 
in which j AB denotes the induced base space metric 

7 AB = 7 ,fl 'gX' ( 58 ) 

whose components will be time dependent, unlike those of the measure ^i ABC . For 
purposes of physical interpretation it will usually be convenient to take this measure 
to represent the number density of ionic nuclei, so that the confined baryon number 
density n c = t u n" will be given by 

n c = A c m , (59) 

where A c is the atomic number, meaning the number of confined baryons (protons plus 
confined neutrons) per nucleus. It is to be observed that, whereas the number current 
will fail to be conserved in the chemical equilibrium case characterised by (|76jl . the 
formalism is such that the corresponding ionic number current 

n{ = m (60) 

must automatically satisfy the conservation law 

V v n? = (61) 

clS db mathematical identity. 



3.2 Chemical gauge adjustments 

For particular physical applications, the specification of the part of the baryon current 
that is to be treated as "confined" will be to some extent dependent on the timescales of 
the processes under consideration as compared with the lifetimes for quantum tunnelling 
through the confining barriers containing the nuclei. For relatively rapid processes the 
fraction of the relevant quantum states that should be considered to be confined may 
be subject to an increase, so the corresponding current, n c u say, of baryons that are 
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effectively confined will be somewhat larger than the value, n v c say, that would be 
appropriate for processes occuring over longer timescales. It is therefore of interest to 
consider the effect of such a chemical basis adjustment, as given by a transformation of 
the form 

n v = (1 + e)< , nt v = < - e< , (62) 

for some suitably specified dimensionless adjustment parameter e. So long as e is just 
a constant, it is immediately apparent that such a transformation will leave the "free" 
4- momentum covector invariant, meaning that we shall have 



7r f „ = tx 1 v , (63) 



and it is also easy to see that in terms of the new variables the stress energy tensor ()56|) 
will also be expressible in the same canonical form as before, 

T\ = T\. (64) 



As in the fluid case jT^J this covariance property of the canonical stress energy 
tensor is not restricted to transformations for which e is constant, but can be seen to 
hold whenever e is specified as a function just of the atomic number A c = n c /nj, and of 
the material position coordinates q A . Since it can be seen from (|57|) that the variation 
of the ionic number will have a "live" (fixed background) part expressible in the form 

Oni 

^ni = th TabCT^^ + q£ 5 1 A , (65) 

it follows that the corresponding variation of the adjustment parameter will have the 
form 

de ( '5^n c . ..... \ ( de de A c dni\ . . . 

= m: (-sr " A ^-^ s < ) + {^-M;7r^) 5q - (66) 

Since we have S^ric = t v 8n" it follows that unlike the free particle momentum, which 
is subject to the chemical invariance condition (|63|h the confined particle momentum 
covector will undergo a non-trivial chemical adjustment given by 

< = < + " (e % + A-^-uH^j , (67) 

while for the coefficient in the final term of (JMj) we shall also obtain a non-trivial 
adjustment given by 

K = K - - * f >c^c^7 AB C^ • (68) 

It can however be seen that the extra terms proportional to de/dA c will cancel out in 
the final expression f|56[) for the canonical stress energy tensor, whose chemical gauge 
invariance property (|64|) is thus confirmed. 
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3.3 Vortex pinning and chemical equilibrium 

The chemical gauge invariance of the "free" (though not the "confined" ) 4-momentum 
covector is essential for the modelisation of superfluidity, whose effect on a mesoscopic 
scale (large compared with the microscopic lattice spacing but small compared with 
intervortex separation) will be expressed by the condition that this covector should 
have the form of a gradient, 7r = KV^ip of a periodic quantum phase angle tp, with the 
corollary that the corresponding chemically invariant (and even Milne invariant [T2"j ) 
"free" vorticity tensor w will also vanish at the mesoscopic level. At the macroscopic 
level, on a scale large compared with the separation between vortices, this vorticity 
will have a large scale average value that will not vanish, but to be compatible with a 
mesoscopic fibration by 2 dimensional vortex lines it must still be algebraically restricted 
by the degeneracy condition 

Kl^U = • ( 69 ) 

This condition will automatically be satisfied by the field equations that will result 
from the relevant variational principle, which will require that the vanishing perturba- 
tion condition (JHJ) should be satisfied, not for independent variations of all four of the 
displacement fields £ f M , £ C M , which in view of the constraint (jH3) would lead to 

overdetermination, but just for variations satisfying the corresponding restraint (|18j). 
This requirement entails the vanishing, not of all four of the relevant force densities / f M , 
/ c , , f s given by the prescriptions (f3~9"j) and (JHJ), but only of the free force density 
f and of the almalgamated force density f 1 given by (fT^j) that acts on the part that 
is convected with the ionic lattice. 

According to the general formula (|39jl . the vanishing of the free force density 
entails not only the separate conservation condition 

V„< = , (70) 

for the free part n f M of the current density but also a dynamic equation of the familiar 
form 

^,< = 0, (71) 
which using (jHHjl can be seen to imply 

n f £tu f = , (72) 

which is interpretable as meaning that the vortex lines are "frozen" into the free fluid 
in the sense of being dragged along with it. It follows, according to (|69|) that there will 
be no dissipation in the sense that the entropy will be conserved, 

V v s" = . (73) 

Instead of postulating the strict application of the variation principle, we can obtain 
a non-dissipative model of another physically useful kind by postulating that the vortices 
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are "pinned" in the sense of being frozen, not into the fluid but into the atomic lattice 
structure, which means replacing (J7T|) by a dynamical equation of the analogous form 

w^u" = , (74) 

which unlike (171)1 is unaffected by the chemical base transformations considered in 
Subsection 13.21 

Whichever of the alternative possibilities (|7T|) or ()74|1 is adopted, the degeneracy 
condition (|69[) will evidently be satisfied, and more particularly the left hand side of the 
identity (|52jl will vanish, which means that in order to satisfy the entropy conservation 
condition (J75j) the model must be such as to satisfy the condition 

(£ c - £ f ) V v n? = . (75) 

In the strictly variational case this requirement will obviously be implemented by the 
separate conservation condition ()70|). which will be physically realistic in the context 
of high frequency oscillations. However in applications to slow long term variations it 
will be more realistic to use a model based on the alternative possibility, namely that 
of chemical equilibrium in the solid's rest frame, as given by the relation 

£ c = 8 l , (76) 

which, like (|74|) can be seen to be unaffected by the chemical base transformations 
considered in Subsection 13.21 

Thus, depending on how we choose between the alternatives (fTTjl or ([74)1. and be- 
tween the alternatives (|7U|1 or (|76|) . we can use the same Lagrangian master function A 
for the specification of 4 different kinds of non-dissipative model, which are categoris- 
able as unpinned with separate conservation or chemical equilibrium, and pinned with 
separate conservation or chemical equilibrium. 



4 Specification of the action 
4.1 External action contribution 

To apply the procedure described above, one needs to prescribe a (primary) equation 
of state specifying the functional dependence of the Lagrangian master function A on 
the relevant independent variables. 

In accordance with the general principles described in the preceding work [TT], the 
relevant Lagrangian will be decomposible in the form 

A = A oxt + A int , (77) 

in which the gauge dependent external part is given in terms of a fixed baryon mass 
parameter m by an expression of the familiar form 

A oxt = -m (n c vl + m v f 2 ) - mn b <p . (78) 
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The confined particle 3-velocity and the superfluid 3-velocity v\ are defined here in 
terms of the corresponding 4-velocities and u\ in the usual way, by setting 



< = n c u v c = n c {e u + v%) , 
and the squares in (fTSj) are defined by 



n» = n f < = m{e v + vf) , (79) 
»f=V«fV. ( 8 °) 



where r]^ is the uniform rank-3 space metric defined for the Galilean frame charac- 
terised by the ether frame vector e v according to the specifications 0- This means 
that the kinetic action contribution is deemed here to be independent a priori of the 
material coordinates q A and their gradients, to which it will however be related "on 
shell" by the application a posteriori of the constraint (J4TJJ) to the effect that the crust 
frame should coincide with that of the confined particles, i.e. 

u v = < , v" = v v c . (81) 

This approach differs from the treatment used for the non-conducting solid limit case 
|llj . in which the constraint (J81|) was imposed in advance, but it leads to results that 
are entirely equivalent on shell. The strategy used here is designed so as to satisfy the 
condition to the effect that 

V^ = 0, (82) 

which has the advantage of ensuring that (as shown in Subsection 14. 5jl it will be fully 
consistent with the formalism that has been developed [T2"l IT3~| ITS'] for the multicon- 
stituent fluid limit case. 

In the formulation used here, the external action provides no contribution at all to 
the extra stress term V^ v so - as in the multiconstituent fluid case ^H] ~ the external 
contribution to the stress energy tensor will be given by the simple formula 

TL» = <P^ + <Pl - m( f> <t» = + . (83) 

in which free and confined kinematic momentum covectors are given respectively by 

Pi = m(r] VIM v^ - \v}t v ) , pi = m(r] Ufl v? - \v^t v ) , (84) 

and the corresponding non-local external momentum covectors are given by 

i^Lv =p i „-m(j)t 1/ , n^ tl/ =p%-m<j>t v . (85) 

Since, according to ()34|) the internal part will not contribute, the 3-momentum den- 
sity (|3*0|) will consist just of the external part, which evidently provides the (chemicaly 
invariant) result 

n v = m(n{Vf + n c v") = mrf^n^ . (86) 
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4.2 Internal action contribution 



It is to be remarked that the Galileian frame independent difference — v v determines 
a corresponding - purely spacelike - relative current vector, 

n£ = rif(Vf — v u ) = rtf — Ufu" , n^t^ = , (87) 

which as well as being unaffected by changes of the Galilean frame is also unaffected 
by chemical base transformations of the form (jfi2j) . which will simply give n±. u = n±. 

The internal contribution A int in (77) has to be independent of the choice of the 
Galilean ether frame vector e M . This means that at a given material base location, 
as specified by the fields q A , this internal contribution will depend only on the scalars 
s = s u t u , rif = n\t v and n c = n v c t v and on the material projections 

n A = = n?q^ , 7 AB = l^t^ > ( 88 ) 

so that its generic variation will be given by 

«,„ = + ^Sn c + d -^6n, + ffSn* + |% 7 « + *^ . (89) 

os dn c oni dn A dj AB dq A 

This provides an associated "connective" variation (in which, with respect to 
appropriately dragged coordinates, both q A and Sq A are held constant) of the form 

, , dA. . , dA. . , dA. . , dA. . , „ dA. . , „„ , , 

4A,„, = -gf 4. + -^f4» c + + 4n r + 3^47- , (90) 

in terms of tensorial coefficients defined by 

dA. t dA. t dA. dA. . 

dnl dn A q ' u ' dY u <9 7 ab J '» q ^ ■ yJ ' 

For the purpose of the present analysis, what we need is the corresponding "live" 
variation, as carried out at a fixed position in a fixed background, for which we obtain 

™» = 9J t Ss+ tf fe <+t? 5 +2 f^ W+|r^ < 92 > 

We now use the constraint u y 8q A u = —q A u 8u u to recombine the terms in such a way as 
to obtain a coefficient "P ^ satisfying the condition ()55|) in an expression of the standard 
form 

dA 

6*A int = -QSs + X c Jn v c + xlK + + ~g^Sq A , (93) 

It can be seen that the required (purely spacelike) value of V ^ will be given by 
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while the corresponding expressions for the internal contributions to the 4-momentum 
covectors will be given by 

<9A. . r OA. . m OA. . f dA. . dA. t 

ds ' <9n c 17 n c dn" ' drtf v dn± 

which implies that we shall have 

It follows, according to the canonical formula (f3l)j) that the internal contribution to 
the stess energy tensor will be given by 

= -0 sH„ + X c v n? + X { X + * K ~ , (97) 

with 

dA OA 
P m = + 2^| 7 - , (98) 

and with the generalised presure scalar $ given, as in the multiconstituent fluid case, 
by 

v& = A int + 6 s - X \n» - X \< . (99) 



The identity (|32jl ensures the symmetric and strictly spacelike nature of the ensuing 
pressure tensor, as given in accordance with (J3T|) by 

p r = T -id av , (ioo) 

which must also be symmetric. It will be obtainable according to the prescription 

p in T = ( s@ - Kx f . + "cXlK - xW) y - , = Y p Y a s pa , (ioi) 

in which the stress contribution is given by an expression of the same form as in 
the non conducting solid case [TT], namely 

dA. , , d(A. Jm) d(A. Jm) , . 

S = 2 — - A. .7™ = 2 m v ' nt/ ; = 2 m ' V ug B „ • 102 

This space projected part can be combined with the associated time projected part, 
namely the comoving energy current defined by 

U?? = -T. *u v , (103) 

int int^ ' V > 

to give back the entire internal stress-energy tensor in the form 

t s = ^r/v - u^tu , (104) 
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in which 7^ is the time dependent covariant metric that is given by 

7^ = 7 A b^<?V ( 105 ) 
where the covariant base space metric 7 AB is such that 

7ab7 BC = £A (106) 

Using (|4*5j) it can be seen that 7^ is rank 3 : 

l„ v u v = 0. (107) 

One can thus interpret this covariant metric as the rank 3 Euclidian metric that is 
obtained by substituing the ether frame vector defined in (jSJ) by the solid's reference 
frame u^. As a consequence 7^ will be such that 

i, P r p = ^-u% = Y,- (108) 

The comoving energy current U Cfl (which would be the same as the ordinary internal 
energy current in a locally comoving frame - meaning a Galilean gauge with e M = 
at the position under consideration) will be expressible in the form 

U c » = u>*U \-?£**n£, (109) 

in terms of the ordinary internal energy U int = —t^T.^ tU e u which can be seen to be given 
simply by 

OA 

U. = — ^ n A -A. .. (110) 

Qua ^ v > 

Note that the first term on the right of this formula is not something that depends 
on the elastic solidity, but something that - except in the absence of the entrainment 
effect discussed in Subsection 14.41 - will still be present even in the fluid limit. This 
means that (contrary to what has been suggested in the literature jHj) in the presence 
of entrainment it will never be permissible to simply identify the internal action density 
with the negative of the internal energy density. 

As in the fluid case, we can go on to construct locally defined material 4-momenta 

pi + xl> ti, = p c u + xl (m) 



.f 



u 



from which, after allowance for the non local effect of gravity, we finally get the total 
momenta for the free and confined particles, 

= rfv -m<j)t v , n c v = \l c v - m(j)t u , (112) 

in terms of which the complete stress energy tensor will be expressible as 

T\ = -0 sH v + 7r f X + TT%n? + W v - V\ . (113) 

with 

V^ = S^ + A. m X + <^. (114) 
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4.3 Elastic energy contribution 

As in the multiconstituent fluid case [TH] it is useful to decompose the internal action 
function in the form 

A int = A ins + A ent , (H5) 
so as to obtain a corresponding decomposition 

^ = S^ + S^ (116) 
in which A , and S u „ are the parts attributable to entrainment, and A. and S u ,, are 

cut en t fj, sr 1 , n s i ns fl 

the static internal contributions that remain when the relative current contributions n A 
are set to zero. For this static part (but not for the rest) the action density will just be 
the opposite of the elastic energy density € as defined by setting 

U ins = €, (117) 

which can be seen from f|110|) to correspond to 

A. =-€. (118) 

ins \ / 



It will commonly be convenient to further decompose this static energy contribution 
in the form 

€ ins = € + € sol , € = -A QJ (119) 

in which €t Q is the part that remains in a relaxed configuration for which € is min- 
imised for given values of the independent current components n c , rif, n A , and of the 
determinant of the induced metric 7 AB . Fixing this determinant is equivalent to fixing 
the value of the conserved number density n\ that is specified by (|57J) . We shall use 
the notation 7 AB and 7 AB respectively for the corresponding relaxed values of 7 AB and 
its inverse 7 AB (as defined by 7 AB 7 BC = 5?) at which, for given values of rij, n c , n f , the 
maximisation occurs. Thus substitution of 7 AB for 7 AB in the solidity term 6* sol or in 
the total € or will give a generically reduced value 

^ ^<^ 5 (120) 

but it will have no effect on the relaxed part, or on the ionic number density rij, for 
which we simply get 

€ n = € Q , ni = m. (121) 



The relaxed contribution will evidently be of ordinary (albeit non-barotropic) perfect 
fluid type with a generic variation of the form 

8€ Q = -5A Q = e Q 5s + x'jnt + X c Jn c + ^Jm - X s QA 5q A , (122) 
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(in which the final term allows for the possibility of built in inhomogeneity in addition 
to the stratification due just to the variation of the atomic number ratio A c ) with 



5n I = in l7AB 5 7 AB + ^^ A . (123) 

According to (f§5|) the relaxed contribution to the pressure tensor (|101|) will be given 
by 

P o B = (°o s + *o nf + ^o nc )^ AB - S o ' ^o" = (^o - X>i)7 AB , (124) 
so one finally obtains an expression of the familiar isotropic form 

P ab = p o7A b ; P =x' n t +x l n I -G , X l Q = X S Q +A cX c Q + -f % • (125) 



Let us now consider the solidity contribution €? aol whose job is to allow for the effect 
of deviations of 7 AB from its relaxed value 7ab (as determined by the scalars n f , n c , m). 
Such deviations can conveniently be accounted for [2] in terms of the constant volume 
shear tensor whose material base space representation is specified as 

?ab = l(7 AB -7 AB ), (126) 

which means that the corresponding space time tensor will be given by 

Sn» = ^abCC = |(T/» - 7/uO • (127) 



In most applications to behaviour of a perfectly elastic (rather than plastic or other 
more complicated) kind, it will be sufficient to use an ansatz of quasi-Hookean type 
[2], meaning one in which the solidity contribution has a homogeneously quadratic 
dependence on the deviation (|126|) in the sense that it will be given by an expression 
of the form 

^ = 5^^. (128) 

with 

^ABCD _ j^(AB)(CD) — PCDAB (129) 

in which j] ABCU [ s the relevant shear elasticity tensor, for which the check symbol is used 
to indicate that, for a given value of the material position coordinates q A , it depends 
only on the scalars s, rif, n c and (via n{) on the relaxed metric 7 AB . The condition 
that 7 AB and 7 ab must have the same determinant entails that on the 3 dimensional 
material base the symmetric shear tensor <^ AB will have only 5 (instead of 6) independent 
components, and more specifically that to first order it will be trace free with respect 
to either the actual metric 7ab or the relaxed 7ab . It is therefore necessary to impose a 
corresponding restriction to completely fix the specification of the solidity tensor 1> BCD ; 
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which can most conveniently be done [2] by requiring that it be trace free with respect 
to the relaxed metric 

i; ABCD 7 CD =0. (130) 

The specification of the solidity (i.e. shear elasticity) tensor is not by itself sufficient 
to complete the specification of the elastic system, as it is also necessary to specify the 
dependence on n\ of the relaxed inverse metric 7 AB . The simplest possibility is that of 
what has been termed a perfect solid |2j, meaning one at which the elastic structure 
at each material position is isotropic with respect to the relaxed metric, which in that 
case can vary only by a conformal factor. This means that it will be given in terms of 
its value 7^ B say at some fixed reference value n say of the ionic number density ni by 

7 ab = (n 7 /n ) 2/3 7 AB , 7 AB = KM) 2/3 7o AB • ( 131 ) 
In the case of a solid stucture that is isotropic (as will typically be the case on a macro- 
scopic scale after averaging over randomly oriented mesoscopic crystals) the rigidity 
tensor in the quasi Hookean ansatz will simply have to be given in terms of the relevant 
scalar shear modulus E by the formula 

f"ABCD _ 2 jp (^A(C^D)B _ I^AB^CD) ^ (132) 

where the scalar E is the rigidity modulus which is usually denoted by /x in the 
litterature, a symbol which is already being used in the general formalism used here to 
design the material momentum components. 

Equation ()128j) will then give the simple formula 

A sol = -€ sol = -E, 2 , (133) 

in which the scalar shear magnitude q is defined by the formula 

q 2 = ^AB 7 CD ?bcW _ I( 7 AB ?ab) 2 ? (134) 

in which the final term will in pratice be negligible since of quartic order, 0{q 4 }, in the 
small <; limit that is relevant, as the trace is already of quadratic order, 7 AB ? AB = 0{<, 2 }. 
Under these condition the solidity contribution will provide a pressure tensor given by 
the formula 

P AB = -£ABOD ?AB + P ol7 AB , (135) 

of which the final term is a pressure contribution given by 

f dE dE dE dE E\ , . . 

which, since it is of quadratic order in will be negligible compared with the first 
(linear order) term for most practical purposes. (It is to be noted that in the final term 
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of (J136j) the sign given here corrects an error in the sign of the corresponding term in 
the relevant equation (6.19) as written in the original treatment [2] of the perfect solid 
model.) The corresponding pressure adjustment contribution for the canonical formula 
(j97j) will be given by 

V* v = 7 „A^ Aw< V " Ui^ + ^) , (137) 
in which, again, the quadratic order term at the end will be negligible in practice. 



4.4 Entrainment contribution 

In general the entrainment action function A cnt will depend, for a given values of q A , on 
the relative current components n A = q A n± as well as on the scalar magnitudes nf , n c , 
and the induced metric components 7 AB , so its generic variation will have the form 

<9A <9A <9A 

5 A = -0 t 6s - y <5n { - y c Sn c + — ^5 7 AB + — ^Sn A + -^Sq A . (138) 

cnt cnt **** 1 ^ c QryAB ' Qua Qq A 



The entrainment action A cnt is characterised by the condition that it vanishes when 
the relative current components n A are set to zero, so when these components are 
sufficiently small, as will typically be the case, it will be a good approximation to take 
this contribution to have the homogeneous quadratic form 

A . = — m X n A n B , (139) 

in which the entrainment mass tensor has components m^ B , that (like the static action 
contribution) are independent of the current components n A , so that the corresponding 
partial derivative in (|138J) will be given by 

dA. t dA t 1 x , , 

_^nt = = _ nB _ 14 Q 

on A on A n.f AB 



It is conceivable that the relaxed action function might involve a built in anisotropy 
favoring relative currents in some particular direction, but in cases of the simplest kind, 
to which the remainder of this subsection and the next will be restricted, this function 
A cnt will be of purely fluid type in the sense that for given values of q A it will depend 
only on the set of five scalar magnitudes consisting of rij, n c , rtf, s , together with the 
relative current magnitude n± that is defined in terms of the (unrelaxed) metric value 
7 AB which - using the material index lowering operation specified by the induced metric 
7 AB - will be given by 

nl = n A n A , n A = 7 AB n B , 7 AB 7 BC = S° . (141) 
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Such a functional dependence provides an expansion 

5A cnt = - Xljm - xL^c - xlM + \L^ A + , (142) 



of similar form to the perfect fluid contribution ()122|) but with an extra term involving 
a partial derivative that provides an expression of the form ()140|) in terms an isotropic 
mass tensor given by 



i 



m AS = < Tab , m\ = 2m , (143) 

while the other partial derivatives in ()138J) will be given by 

9A nt i fm l c , s \ <9A cnt s s <9m 



djAB 2 ^ „ f "'A^B + Knt n ll^j > ^7 A entA Xnt £ gA ' ( 144 ) 

The scalar m f c introduced in this way is identifiable as the increment m f c = m* — m of 
the effective mass of the free baryons (meaning the superfluid neutrons) as compared 
with the ordinary baryonic mass m. This mass increment is expected to be positive 
(and in some layers large) ^Bj in the solid neutron star crust, but (moderately) negative 
in the fluid layers below. 

4.5 Relaxed action contribution 

It will be useful for many purposes to replace the decomposition (jll5j) of the internal 
action density by an alternative decomposition of the form 

A int = A lax + A sol , (145) 

in which the relaxed - meaning shear independent - part will evidently consist of the 
combination 

A lax = A G + A cnt • (146) 

The use of such a combination is particularly convenient whenever the entrainment 
contribution is of the isotropic type characterised by the variation expansion ()142l) . 
in which case the complete relaxed action density will have a variation given by an 
expansion of the analogous form 

f 

777 

5A = -0 Ss - Y 5n { - X c Sn c - X S Sm + X s 5q A + —^Snl, (147) 

lax lax A, lax I A, lax C A, ent 1 lax A ^ 

with 

9 =0 +0 f y f = Y f + Y f y c = y c + y c 

W lax ' cnt ' A lax A Acn t ' A lax A Acn t ' 

XL = X S +Xt, Af axA = A s OA + A c S ntA . (148) 
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It can be seen that the relaxed contribution to the "live" action variation (|93|) will 
simplify to provide an expression of the form 

c^A = -0 5s + v f v 5n? + x c v S< ~ X S ^n T + A, s 4 6q A , (149) 

lax lax 'Max 1 ' 1 'Max" c ^ lax A - 1 " \ ' 

in which v f „ and y c ,. can be read out as the relaxed parts of the internal momenta 

' Hax ^ 'Max 1 ' A 

given by J55J). Since we can write S^rij = t^S^rij 1 , it can be seen that, as the analogue 
of these free and confined particle 4-momentum contributions, we shall also be able to 
read out a corresponding ionic stratification 4-momentum contribution x S v i so as to 
obtain a complete set of relaxed internal momentum covectors that will be given by 

f f 

TfL Tfl 

xL = -xL^+ ^^ nl ' xL* = ^^ nl ' xL» = -*J» ■ ( 15 °) 

In the formula (|9*4^) for the extra stress, it transpires that the contributions involving the 
mass increment m { c (proportional to <9A lax / 'dnl) will cancel out, leaving the expression 

K„ = -xL n n AB q^, (151) 

in which the part due to stratification is all that remains. It can be seen that the 
corresponding space time tensor will be given simply by 

^ = - *V ■ ( 152 ) 

In terms of their contributions to the total momenta, which will simply be given by 

^laxf = _ Qax^f > ^lax" = XlaxV 1 71 lax" = Xl^V 1 71 lax" = X lax" ' (153) 

the corresponding elastically relaxed force contributions acting on the entropy current 
and on the free and confined particle currents will be given by an ansatz of the standard 
form (|39~j) which gives 

t = 2^] + iV M ^, (154) 

C = 2 <V[^U + < X - V X > ( 155 ) 

and 

/I, = 2 <V [M <^ + n^VX , (156) 

while for the analogously defined force contribution (J4*U|) due to stratification acting on 
the underlying ionic lattice, it can be seen that similar reasonning leads to an expression 
of the slightly different form 

fL = + \i u , \i u = ■ (157) 
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In terms of the amalgamated ionic 4-momentum contribution defined by 

= + A c <^ + (s/m)^ , (158) 
the corresponding contribution 

to the amalgamated ionic force density (fT9|) can be seen to be expressible directly by 
the formula 

fL = 2 <V [M 7T 1 I axi/] + < tt^V^ - m^V, ( a /m) + \ s x . , (160) 

while the associated contribution to the stress energy tensor (J3f))) will be given neatly 
by 

T " = (A, - n/V fl - n/V 1 J 5^ + n/V „ + n/V „ . (161) 

lax^ V lax t laxP 1 laxP' " I lax^ 1 lax^ V > 



4.6 Complete description for perfect conducting solid. 

Replacing the relaxed contribution ()158|) by the corresponding total ionic 4-momentum 
covector given (in terms of the confined atomic number A c ) by 

ti\ = nl + A c ir c u - (s/m)0 t v , (162) 

we can apply an analogous tidying up operation to the complete stress energy tensor 
(|36p which will thereby aquire the form 

T\ = (A - n f V p - n{n l p )5^ + + - V* v , (163) 

in which the only manifest allowance for the effects of solid rigidity is in the final term. 
In the case of a perfect conducting solid - meaning one whose structure is fully isotropic 
with respect to the relaxed metric - this final term v will be given by the formula 
(I137J1 which will be expressible to first order in the shear amplitude q by a prescription 
of the simple form 

V» lV = 2E 1 ^ pv + 0{e}, (164) 

in which E is the ordinary shear modulus (which elsewhere is commonly denoted by 
the symbol \x that, in the present context, has already been used for the designation of 
the material momentum components). In principle, the other terms in ()163|) will also 
be influenced by the solid rigidity contribution A sol , but since such effects too will be of 
quadratic order in they will be effectively negligible in applications of the usual kind, 
in which deviations from an elastically relaxed configuration are small. This means 
that for practical purposes, in a conducting solid of the perfect (meaning intrinsically 
isotropic) type, the deviation from behaviour of (multiconstituent) fluid type [12] will 
be entirely contained in the extra term given by (|164j) . 
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5 Non conservative generalisation 



5.1 Dissipative interpolation 

In the non-dissipative models described above, it can be seen that the introduction of 
the entropy density s as an indpendent variable was in practice redundant, since its 
effects could be allowed for simply by a readjustment of the stratification, because the 
ratio s/rii was fixed on each material world line so that it depended only on the material 
base space variables q A 

The reason for taking the trouble of introducing it is that the entropy density will 
acquire a non-trivial role as soon as these convectively conducting solid models are gen- 
eralised to allow for dissipative effects in the manner that has recently been described 
in detail ^H] for the fluid case. The second law of thermodynamics tells us that, in a 
dissipative application, the entropy current s v need not satisfy the conservation condi- 
tion (J75|) . but that - for a system that is closed in the sense of being thermodynamically 
isolated from the rest of the universe - the model must be such as to ensure satisfaction 
of the inequality 

V^>0. (165) 

On the assumption that it is isolated not just thermodynamically but in the stronger 
mechanical sense of (146)1 . we have seen that a system characterised by an action of the 
kind presented above, and thus by a stress energy tensor of the form ()56|) . must automat- 
ically satisfy the energy conservation identity (J52|) . Subject to the usual presumption 
that the temperature is positive, this means that the second law requirement (|165|) 
will be expressible for such a model as 

(S c - £ l )V u n? + n/V^X > • (166) 

The simplest and most obviously natural way of satisfying this positivity requirement 
starts by postulating that the two terms in (|166|) are separately positive. For the first 
term this leads to ansatz to the effect that relevant neutron fluid creation rate should 
be given by an expression of the form 

V„< = k(£ c - S { ) . (167) 

in terms of some positive transfusion coefficient k that might be expected to be a 
sensitive function of the temperature, 0, on which the weak interactions that would be 
involved are known ^Hj to be highly dependent. (It is to be noticed that the formula 
(1167)1 will be preserved by the chemical gauge transformation (JH2*|) for any fixed value 
of the adjustment parameter e.) 

The treatment of the second term in (|166)1 is not so simple, because it is necessary to 
respect the degeneracy requirement (|6T?j) that ensures that the mesoscopically averaged 
vorticity w i is orthogonal to two-surfaces that represent the flux of quantised vorticity 
tubes with generators of the form u[ + V[ for some spacelike vector V[ such that 

+ V?) = , V?t„ = . (168) 
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In terms of such a vector, the positivity condition on the second term in will be 

expressable (assuming positivity of ri{) as 

u^ v Vr>0, (169) 
so a vector of the required form will be given by an expression of the form 

^ = -^ f /, (170) 

in which c r is a positive drag coefficient that has been adjusted so as to be dimen- 
sionless by the inclusion of the denominator «/ f , which is defined as the magnitude of 
the (spacelike) vorticity vector u> f/i = \e >iVp 'uj i vp . This provides a superfluid equation of 
motion of the form 

u?m l pal ^ = c r A M , (171) 

in which v r ^ is a relative flow velocity, namely that of the medium relative to the 
vortex sheets, as defined in terms of the relevant orthogonal projection operator by 

v T " = , -It = (w* )" W^l^V • (172) 

On the basis of work by Jones fH] , a formula giving a rather low value for the required 
drag coefficient c r , in the low temperature limit, as a function of the relevant densities, 
ni, n c , rif, has been provided by Langlois et al. but it depends on microscopic 
parameters that are difficult to evaluate, and its validity is in any case a subject of 
controversy. In the very different picture developped by Alpar et al. (201 ^ has been 
suggested that the result will be highly temperature dependent, and that the effect 
will be describable as (nonlinear) "creep" rather than ordinary (linear) drag, in the 
sense that the appropriate coefficient c r will not just depend on the relevant scalar 
densities as well as the temperature 0, but that it will also be strongly dependent on 
the magnitude v T of the relative velocity ()172|) (according to a formula of the form 
c r oc v T _1 arcsinh{f r /i) r } for some velocity independent - but temperature sensitive - 
quantity v T ). 

It is to be remarked that the chemical reaction rate formula ()167j) provides an in- 
terpolation between the non dissipative limits k — > oo, namely the thermal equilibrium 
case, and the variational case n — > 0, for which free and confined particle currents are 
separately conserved. Similarly the drag or creep formula (|171|) provides an interpola- 
tion between the non dissipative limits c r — * oo, namely the perfect pinning case, and 
variational case c r — > 0, for which the vortices are freely transported by the superfluid 
flow. 
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5.2 Conclusion 



Subject to the prescription of two suitable secondary equations of state for the coef- 
ficients k and c r as functions of the relevant variables, the foregoing equations (J167|) 
and (|171|) constitute a complete system of equations of motion for the nine indepen- 
dent component variables (which can be considered to be the four superfluid current 
components n u , and the five materially convected components q A n c and s) when used 
in conjunction with the baryon conservation condition (j51|) and the condition (j4T?j) of 
conservation of the stress energy tensor T^ v that is obtained, according to (j!63J) . from 
the the master function A, whose prescription, by a suitable primary equation of state 
is described in Section EJ 

The stress energy conservation condition (j4Tlj) expresses the requirement that the 
system should be effectively isolated, not just mechanically, but also thermally. This 
last requirement will not be entirely realistic when the transfusive adjustments governed 
by (j!67)) are taking place ^Hj, since the beta processes involved will create neutrinos for 
which the stellar medium will be effectively tranparent so that instead of being locally 
confined they will rapidly escape from the system. The ensuing heat loss can easily be 
formally taken into account ^3] by replacing ()46|) by a generalisation in which there is 
an extra term so, that it takes the form 

V ll T» v = -pV v <l>+Qt ll , (173) 

in which Q represents the heat loss rate per unit volume. This means that the energy 
conservation identity (p)2*|) will need to be replaced by the energy loss formula 

Q + = (£ c - £ f ) V„< + nfzu^X = « (£ c - £ f f + c r Wf w ! v T 2 , (174) 

in which the right hand side is the combination of terms that (by the non negativity 
of the coefficients k and c r ) has been made to satisfy the positivity condition Ijlfifij) . 

If Q is prescribed (by what would be a third secondary equation of state) as a 
function of the relevant densities - particularly that of the entropy which determines 
the temperature G - then (J174|) will in principle provide what is needed for calculating 
the evolution of the entropy current. However, as the temperature dependence of Q is 
likely to be highly sensitive ^H], such an approach might not be accurate in practice. It 
might be more realistic to suppose that the temperature would adjust itself to the rather 
low roughly constant value needed to avoid accumulation of entropy, which incidentally 
is the consideration that justifies, as a reasonable approximation, the neglect of heat 
conduction in this work. This means adopting the conservation postulate (J73*Jl to the 
effect that the second term on the left of ()174j) will simply drop out. The implication 
is that Q will be given just by the positive terms on the right of ()174j) . and that it will 
therefore vanish in the non dissipative limit cases that are obtainable, as described in 
Subsection 13.31 by taking infinite or zero values of the dissipation coefficients k and c r . 
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